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Abstract 

The thermodynamics of a gas of strings and D-branes near the Hagedorn transition is 
described by a coupled set of Boltzmann equations for weakly interacting open and closed 
long strings. The resulting distributions are dominated by the open string sector, indicating 
that D-branes grow to fill space at high temperature. 



1 Introduction 



Strings exhibit interesting behavior at high temperature. A gas of closed fundamental strings 
approaching the Hagedorn temperature Th enters a long string phase, where any added 
energy goes into forming string rather than raising the temperature beyond Th- In the 
long string phase the average number of long strings of a given length i is proportional to 
l/£ and the average total number of long strings is log{E), where E is the total energy 
in the microcanonical ensemble. These results are readily found by solving a Boltzmann 
equation for long closed strings |]l], § . They can also be obtained from the exponentially rising 
density of states for free strings, provided one is careful to properly define the microcanonical 
ensemble [||, 0, 3. At very high energy density the description in terms of weakly coupled 
strings in flat spacetime breaks down and the physics of the system is unknown. 

Strings are not the only extended objects in string theory and it is natural to ask how the 
above physical picture is changed when the system also contains p-branes. We will address 
this question in the context of Type II strings and D-branes, by extending the Boltzmann 
equation approach of ref. to include the auxiliary open string sector that describes the 
D-brane dynamics. We find that the long closed string phase is suppressed in the presence 
of a dilute gas of D-branes (and anti-branes) and the system is instead dominated by long 
open strings as the Hagedorn temperature is approached. 

The interactions of fundamental strings always include gravity and the thermodynamic 
limit, where one considers a fixed energy density as the volume becomes arbitrarily large, is 
not well defined for gravitational systems due to the Jeans instability. For a fixed value of 
the string coupling g we can only consider a thermal ensemble on finite length scales. 



where p is the energy density. Alternatively, the limit of large volume can be taken for fixed 
p only if g is simultaneously scaled to zero fast enough. 

This has important consequences if we wish to include D-branes in the system. A Dirichlet 
p-brane carries a conserved charge that couples to a Ramond-Ramond, rank p+1, antisym- 
metric gauge potential. A finite volume target space must carry zero net charge, and as 
the D-branes are the only objects in the theory that carry the R-R charge, we must have 
an equal number of branes and anti-branes in our system.Q The annihilation of branes and 
anti-branes will deplete their numbers in a true equilibrium configuration unless the energy 
density is sufficiently high, p ~ l/gf, to have brane anti-brane pairs forming at a rapid rate. 
At such high energy densities, however, the dynamics of D-branes are unknown and our 
approximations break down. 

We will instead consider a dilute collection of branes and anti-branes interacting with 
a gas of strings. This is not an equilibrium configuration but if the average time between 

"'^Some of the p-branes may be in bound states with (p+2)-branes or (p+4)-branes but the net p-brane 
charge still has to be zero for all p. 
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brane anti-brane encounters is long compared to the time it takes the string gas to reach 
equihbrium, then the configurations at intermediate time scales will be well approximated 
by solutions to our transport equations. Note that we cannot take the D-branes to be too 
dilute. At any fixed value of the string coupling, the Jeans instability limits the size of 
the entire system and we want the average separation between branes to be small on that 
scale. As we shall see below, for weakly coupled strings there is a range of D-brane number 
densities where both of these requirements are satisfied. 



2 Boltzmann equation for closed strings 

We use Boltzmann transport equations to obtain the distribution of string lengths, both 
for the closed strings in the bulk of spacetime and the open string sector associated to the 
D-branes. Boltzmann equations for strings have previously been studied in a discrete model 
by Salomonson and Skagerstam [0, and in a continuum approach by Lowe and Thorlacius 
which can be easily adapted to account for the two coupled string sectors. 

To review, let us first consider the simpler system of only closed strings and no D-branes. 
The Boltzmann equation for weakly coupled closed strings is 0, 

^ = ^{-l<^'f{i) + y^di' ni-i')fii')fii-i')-ifii)fdi'i'fii') 

+ rdi'ii+nfii+n}, (2) 





where f{i) is the average number of strings of length £, k is a positive constant proportional 
to g'^, and V is the volume of the system.^ The terms on the right hand side all involve the 
three closed string interaction shown in Figure la. Higher order interactions of four or more 
closed strings may be ignored at weak string coupling. 

The first term corresponds to a loop of length i self-intersecting and splitting into two 
strings of length £' and i — while the second term describes the reverse process where two 
smaller strings join to form a string of length i. The third term is due to strings of length 
i and i' joining to form a single string of length i + i', while the fourth term represents a 
string of length i + i' splitting into strings of length i and i'. 

In each of the four terms, two string segments meet and exchange ends at the intersection 
point. The interaction rate is proportional to the coefficient k which involves a sum over 
relative momentum and orientation of the segments. In the long string limit this sum does 
not depend on string length nor on whether the two segments belong to separate strings or 
the same one, and thus the coefficient is the same for all four terms. The overall factor of 
1/V refiects the fact that the string segments have to be at the same place in the embedding 
space in order for an interaction to take place. 

At equilibrium we are interested in a static solution of the Boltzmann equation, where 
§lf{() = 0, and the terms on the right hand side of (g) cancel. The equilibrium solution is 

^We have set a' = 1 for convenience. 
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Figure 1: Leading order string interactions. 



easily obtained 0, 

= -f'"' , (3) 

where L is the average total length of string in the ensemble. 

The energy of a long string is to a good approximation proportional to its length, e = ai, 
where a is the string tension {a is 0(1) in our units). The distribution /(£) is related to the 
single-string density of states uj{e) as follows: 

f{i)di = uj{e)e-^'d6. (4) 

Inserting the equilibrium solution (j^) we read off 

u{e) = (5) 

where 

(3H = f3--, (6) 

is identified as the inverse Hagedorn temperature. The corresponding multi-string density 
of states is 

n{E)=exp{l3HE), (7) 

and the single-string distribution function d{e, E), which gives the average number of strings 
carrying energy between e and e + de in a. system of total energy E, is given by 

This result is valid foTC<6<E — c where c is some constant independent of 0]. 

The average total number of long strings is then log-E. At an energy density of order 
one in string units the average total length of string in the gas is L ~ so the typical long 
string has length 

V 

(9) 



logV^ ' 

which is large compared to the linear size of the system, Lq = V^l'^ 
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3 Boltzmann equations for strings and D-branes 



We now include a dilute gas of D-branes in the system. We consider parallel p-branes and 
anti-p-branes in a target space = R®/C|| (8)/C±, where /Cy and /C^ have dimension d\\ = p 
and d± = 9 — p respectively. For simplicity, we will assume < p < 6. The branes are 
wrapped around /Cy, which we will take to be compact with volume of order one in string 
units, but are free to move around in the transverse space }C±. The transverse volume V± is 
macroscopic but bounded by the Jeans volume, 

v'x<y . (10) 

Let n be the number density of p-branes and anti-branes in the transverse space, 

N, + N, _2N, 

The typical separation between branes is Ip ~ n~^/^-^. We must take this distance to be large 
in string units for otherwise the energy density due to D-branes would no longer be small 
compared to 1/g. 

The string gas now contains both closed and open strings. Let /(£) denote the average 
number of closed strings of length i in the system and p{i) the average number of open 
strings of length £. There are two coupled Boltzmann equations which express the rate of 
change of f{i) and p{i) with time. For weakly interacting strings the terms in the Boltzmann 
equations correspond to the various string interactions shown in Fig. |l]. Let us consider the 
f{£) equation first, 

/■°° -I b 

+ di'ii + e)f{i + f )} + ^P{i) - anifii) 

+-{ di\i' - e)p{i') - if{e) di'i'pin}- (12) 

The first four terms on the right hand side correspond to the closed string interaction in 
Figure la and are the same as in the Boltzmann equation of the previous section. The 
remaining terms correspond to closed strings interacting with open strings attached to D- 
branes, as indicated in Figure lb and Ic. A splitting interaction where a closed string is 
converted into an open string, or an open string splits into two open strings, comes with 
a factor a, which is proportional to g and includes a sum over the momentum and relative 
orientation of the string segments involved. The reverse joining process is characterized by 
another coefficient b, which is also proportional to g. 

The equation for p{i) receives contributions from the open string interactions in Figure Id 
and le, and also from the mixed interactions in Figure lb and Ic, 

^^^^^ - ^ I'' d£'p{e')p{£ - t) - anipie) + 2an H de'p{i') - -^p{i) H de'p{£') 



dt 2nV Jo Je nV 
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+^j diMh)J di2p{i2) - ip{i) J di'i'p{i')} 

dn'f{i')} + anif{i)-^p{i). (13) 
At equilibrium, the distribution of strings does not change with time, 

The coupled equations (|1^) and (0) appear considerably more comphcated than the 
Boltzmann equation (^ for closed strings. One nevertheless finds a simple equilibrium 
solution, 

f{£) = iexp[-£/L,], p{£) = ^exp[-£/L,]. (15) 

The distributions are characterized by the average total length of closed and open string and 
the average total number of open strings, 

roc noc roc 

Lc= / diifii), Lo= / d^^p{^), No= / dip{i). (le) 

JO Jo Jo 

These parameters satisfy the relations 

No 2an^V ,r r r 

Lc b 

It immediately follows that the open strings tend to dominate the ensemble. Whenever there 
is a macroscopic number of open strings present. A*";, >> 1, the total length of closed string 
will be vanishing compared to the total length of open string, Lc « Lo- In other words, the 
D-branes efficiently chop up the strings so that the closed string gas is replaced by a web of 
open strings attached to the D-branes. 

The typical length of a long open string can be obtained from the single string distribution 
function for open strings do{e, E) in the microcanonical ensemble. The first step is to read 
off the single-string density of states for open strings from the equilibrium solution ([T5|) , 

ujo{e) = ^"[^ ^ exp (Phs)- (18) 

OCT 

The corresponding multi-string density of states can be obtained in a saddle point approxi- 
mation, which is valid when the total number of D-branes in the gas is large. 



noiE)^e.v{f-^ + PHE}. (19) 
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The open string distribution function is then 



do{e,E) 




ha 

The distribution function has a characteristic energy scale 



(^-f^)\ (20) 




(21) 

below which the distribution is more or less fiat and above which it is strongly suppressed. 
The energy of a long string is proportional to its length so there is a corresponding charac- 
teristic length scale for the open strings in the gas. Since the rate coefficients a and h are 
both proportional to g the characteristic length is given in string units by 

/o ~ (22) 

n 

up to factors of order one. At the Hagedorn energy density, p ~ 1, the characteristic open 
string length is long compared to the typical D-brane separation Ip ~ n~^^^^, so the D-branes 
are effectively overlapping. 

Our calculations have involved a number of assumptions and we close this section with a 
discussion of their range of validity. The Boltzmann equations only have the relatively simple 
form of ([121) and (|13D if we can assume that points on a string that are separated by a finite 
parameter distance are at uncorrelated positions in the embedding space. This criterion is 
satisfied if the volume Vg occupied by a long string when placed in an infinite space is large 
compared to the system volume V. The embedding of a long string of length £ is a random 
walk which occupies a volume Vg ~ i'^^/'^ in d± spatial dimensions. Let us assume that the 
energy density in the system is order one in string units and that the spatial volume is close 
to the maximum value allowed before the Jeans instability leads to gravitational collapse, 
V ~ g~'^-^. The characteristic open string length (p2D involves the D-brane number density 
and we find that Vg > V is satisfied only if 

n < g\ (23) 

We also want to impose the condition that the typical D-brane separation, Ip ~ n"^/*^^, is 
small compared to the linear size of the system, which requires n > g'^^. Since (i_L > 3 for 
p-branes with p < 6 there is always an allowed range for the number density. 

If Vg for typical strings is smaller than the system volume the terms in the Boltzmann 
equations that correspond to string self-interactions have to be modified by replacing the 
factor of V in the denominator by Vg. We have not found explicit solutions of the resulting 
equations but we expect such a modification to increase the weight of shorter strings in the 
distribution when the D-brane number density is higher than in (^3]). 
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Our results also rest on the assumption that the strings take less time to reach an equi- 
librium configuration than the timescale on which the D-branes and anti-branes meet and 
annihilate. It turns out that this is a weaker condition than (0). The D-brane lifetime can 
be estimated as 

r , (24) 

nav 

where a is the brane anti-brane annihilation cross section and v is the average D-brane 
speed. The D-branes are immersed in a string gas with an energy density of order one in 
string units. Their kinetic energy will therefore be of order one which means that their 
motion is non-relativistic with v ~ g^^'^- The eikonal approximation to the low velocity 
scattering amplitude of a D-brane and anti-brane pair develops a tachyonic instability if the 
impact parameter is less than a certain value of order one in string units p . We use this as 
an estimate for the annihilation cross section and take a ~ 0(1). The D-brane lifetime is 
then found to be r ~ g^^^'^n'^. 

This is to be compared to the time it takes the string distribution to find equilibrium. For 
an ensemble dominated by open strings the relevant processes are the open string interactions 
in Figure Id and le and we can ignore the closed strings. In the long string limit the 
probability for a pair of string segments, that happen to come across each other, to interact 
is independent of the overall string lengths. Furthermore, if a long string undergoes an 
interaction its new length will be largely uncorrelated to the old one. It follows that a gas 
of long strings will find equilibrium on a time scale where the strings have interacted a 
few times each. It turns out that the exchange interaction Figure le, which occurs in the 
bulk of spacetime, is more efficient at redistributing length among the open strings than the 
endpoint interaction in Figure Id, which can only take place inside a D-brane worldvolume. 
The probability for a particular string to participate in such an exchange interaction in a 
unit of time goes as g'^i, where i is the length of the string. Using the characteristic length 
(P^), and assuming p ~ 1, the time to reach equilibrium is 

re - p. (25) 

It is short compared to the D-brane lifetime if the D-brane number density satisfies n < g^^'^, 
which is a weaker condition than (^). The equilibrium time scale Tg can also be estimated 
by linearizing the long string Boltzmann equation (|13D around the equilibrium solution (|15D 
and find the rate at which the perturbation decays. This second method gives an answer 
that agrees with (p5|). 



4 Discussion 

D-branes strongly affect the long string phase associated with the Hagedorn transition. In 
the absence of D-branes it consists of long closed strings that traverse the entire system 
many times over, but a dilute collection of D-branes will suppress the long closed strings and 
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replace them with long open strings that are attached to the D-branes. The typical open 
string is long compared to the D-brane separation and is therefore most likely to have its ends 
on different D-branes. Since open strings are associated to the dynamics of the D-branes 
themselves these results suggest a physical picture where D-branes grow with temperature 
as the Hagedorn transition is approached and eventually overlap to fill space. 

The system enters the long string phase when the energy density becomes of order one 
in string units. For a D-brane number density that satisfies the condition the 'bare' 
energy associated with the D-branes is then small compared to the energy carried in the 
open strings attached to them. If the energy density of the system is increased further the 
characteristic open string length increases and the D-branes become more interconnected. 
If the energy density becomes of order 1/ g non-perturbative effects, such as pair creation of 
branes and anti-branes, become important and our description in terms of weakly interacting 
strings breaks down. 

In this paper we took the string coupling to be weak and limited the system size in 
order to avoid the Jeans instability but it would be interesting to incorporate long range 
gravitational effects into the description of the string distributions. Some related issues are 
addressed in recent work of Horowitz and Polchinski . 
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